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Abstract 

In the article by M. G. Mustafa published in Phys. Rev. C 72, 014905 (2005) the author has estimated 
the total energy loss of a charm quark and quenching of hadron spectra due to the collisional energy loss of 
energetic partons in an expanding quark-gluon plasma employing Fokker-Planck equation. We wish to point 
out through this Comment that some of conceptual and numerical results of the said paper are unreliable. 
For the sake of clarity our discussion will focus on the massless case (although a few remarks on the m / 
case are also made). 
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General remarks 



The said paper by Mustafa and an earlier article by Mustafa and Thoma [0] were aimed at 
expanding further the pioneering works by Baier et al., yl] and Miiller [0] on the important phe- 
nomenon of jet quenching, i.e., the suppression of hadronic spectra produced by partons (having 
large transverse momentum) which suffer sizable energy loss while traveling in a hot and dense 
quark-gluon plasma. Refs. [[l], 0] attempted to fulfil this aim by explaining the relevant physical 
mechanisms, recalling the expressions of coUisional & radiative energy loss rates predicted by 
finite-temperature QCD, setting up a Fokker-Planck transport equation for massless or massive 
partons, solving the partial differential equation analytically by using Fourier transform together 
with the method of characteristics, and plotting the predicted quenching factor Q{pt) and Z^/vr ra- 
tio graphically for a Bjorken expanding QGP expected to have been produced in ultra-relativistic 
Au+Au collisions at RHIC. Unfortunately, as we show below through the present Comment, there 
are many questions which are either ignored or not clarified properly in the works of 

Notations 

Although the emphasis of [1] was on the massive charmed quark yet our discussion of vari- 
ous controversial questions becomes easiest in the massless case. Employing c = 1 units let the 
leading parton moving in the transverse, i.e., x direction be associated with mass m = 0, initial 
momentum pq, initial energy Eq =\ pq \, travel time t, distance covered L = t, residual momen- 
tum variable p = po — 5p, residual energy variable E =\ p \= Eq — e, coUisional momentum loss 
variable 6p, energy loss variable e, drag coefficient A(t), diffusion coefficient Dp{t), useful inte- 
gral G{t), average momentum {p)t, useful kinematic symbol r]{p, t), and time dependent variance 
Wit) defined by 

G(t) = ea:p{-/o*rft'A(t')}, r,=p-{p)t [CI) 
W{t) =AG^{t) /o* dt' DF{t')/G^{t'). (C2) 

The first lacuna of ^ 0/ is that they contain mix-up of notations like p + e at several spots [e.g., 
eq. (3) of ref. (1)] and also they fail to mention the fact that E ^\ p \, e ^ 5p, L ^ t in 
the massive (m ^ parton case. One may argue that, since jet quenching is a high-energy 
phenomenon with p > 5 GeV typically, hence the charm quark mass m ^ 1.5 GeV can be 
ignored. This argument, however, is untenable because the fundamental quantity derived in Jit] 
was the probability distribution Di{p, t) which must be computed even down to p = and plotted 
from p = 1 GeV onwards [Fig. 4 of ref. 1]. In these regions the charm quark mass cannot be 
neglected, i.e., the care E ^ \p\ and L should have been exercised in yj]. 

Fokker-Planck equation 

Refs. IE 0] are based on the fundamental transport equation [see eq. (35) of ref. 1] 

^ = A|(pZ},) + D^^ ; 0<t<oo, (C3) 

where the transport coefficients are treated as momentum-independent according to the classical 
view and a suffix / has been attached to the one-body distribution -D(p, t) for the sake of clarity. 



Imposing 6{p — po) initial condition the unique solution of (C3) is found to be [see eq. (42) of ref. 
1] 

= ' -^<^'<^' (C4) 

where the exponential contains purely the second power of r]. Unfortunately, the conventional 
version (C3) of the FP equation cannot properly describe relativistic brownian particle. The (1 + 1)- 
dimensional covariant Langevin dynamics [5] leads to a new FP equation in the massive case as 

where Dnip^t) is the Hanggi-Dunkel distribution function, A(t) and B(t) may be taken as 
momentum-independent coefficients, and the replacement E — » \p\ — > 1 77 1 is to be done in the 
massless limit. The second lacuna o/iH 12/ is that they did not employ the relativistically logical 
transport equation (C5) and also did not seek different type of time-dependent solution Dh in 
which the exponential may contain \ t] \ linearly. 



Probability function in p versus E 

Even if the solution (C4) in momentum variable is assumed to be valid yet care must be ex- 
ercised while transforming to the energy language. We know that both the points ±p are mapped 
into the same E =\ p \. Hence, by standard mathematical statistics, the probability distribution in 
E would read 

Dii{E,L) = Di{p,L) + Di{-p,L) ; < E < 00. (C6) 

which is automatically even in p and preserves the total probability in the sense that 
JZo dpDjip,t) = C dEDniE,L) = 1. 

Unfortunately, for non-zero {p)t, the piece Dj{p,t) = Djjj(E,L) itself is not even in p, it 
cannot represent the probability function over the positive E axis, and its normalization integral is 
awkward namely 

J,^ dpDj{p, t) = J,^ dE Djjj{E, L) = i erfc (^) {C7) 

with erfc being the complimentary error function [6]. The third lacuna of^ 12/ is that they wrongly 
assumed the energy-space probability function to be simply Djjj{E, L) = Di{p, t) [see eq. (28) 
of ref . 2] which contradicts the logical expression (C6) and also yields awkward normalization 
(C7). 



Unphysical asymptotic form 

Suppose the philosophy of is applied to a thought situation where the massless Brownian 
particle is embedded in a huge thermal bath at finite ambient temperature T^o ■ Keeping in mind 
the definitions (C2,C4) we expect that at instants large compared to relaxation time 

G{t)^0 ■ {p)t^O ■ W{t)^W^, iC8) 

^Kp, t)-^7^exp{^}^j^exp{^}, (09) 

The fourth lacuna 0/ /til 0/ is that the equilibrium distribution of their test particles does not tend 
towards the massless Maxwellian form as t ^ 00 even if Bose/Fermi statistics are ignored. It is 
known that the correct asymptotic Maxwellian is retrieved if stationary solution is sought for the 
Diinkel-Hanggi js]] version of the FP equation (C5). 



Mean energy and loss 



Suppose the energy space distribution function for < -E < cxd is just assumed to be 
Dni{E, L) = Dj{p, t) (which it is NOT). Then, by making the transformation 



(T = {p-{p)t)/VW ; ao = -{p)t/VW, (CIO) 

the average energy of the Brownian particle would read 

{E)ni = r dppDj{p,t)/J,^ dp DM = 
{{p)t erfc{ao) + ^W/n exp{-a^)] /er/c(ao), 

where erfc stands for the complementary error function [@]. Clearly, {E)iii ~ {p)t ~ PoG{t) 
only if I (To | = | {p)t/VW |^ 1. The fifth lacuna o//til 0/ is that, instead of( Cll), they report the 
average energy as simply Eq G{t) [see eq. (31) of ref. 2] without mentioning any restriction on 
the ratio | {p)t/VW \. Of course, if the correct distribution Djj{E, L) defined by (C6) were used 
the corresponding mean energy would become {E)ii = Eq G{t) rigorously. 



Convolution integral for hadron spectrum 

Let the shorthand notation T.{E, L) = dN/cPpx denote the Pt distribution of hadrons pro- 
duced by a jet of energy E after covering a distance L. Then, assuming again the probability 
distribution to be Djjj(E, L) (which it is NOT) the medium modification over the vacuum spec- 
trum is given by 

T.^-\E,L) = /o°° deeDui{E,L,Eo = E + e) S--(Eo). (C12) 

which is a function of two variables E and L since e has been integrated out. Employing a linear 
Taylor expansion of T.'""-^(E + e) around E Mustafa-Thoma [see eq. (3) of ref. 1] arrived at the 
estimate 

E^^^ f« S'''^^(E + AE) ; = ^o{l -G'(L)}. (C13) 

We wish to point out that the above estimate is unreliable because the experimentally 
parametrized light hadron production cross-section at RHIC energy namely iQ] 

S'''^^TOcx(l + ^)-\ (CU) 

is a rapidly varying function of Eq. It is much safer to compute the integral (CI 2) by Simpson 
quadrature which reveals that the ratio R = Ti^f^^ /Tisimpson^ for E = 5 GeV, is close to 1 at 
L = 1 fm but becomes about 170 at L = 10 fm. (see Fig. 1 for a more detailed view). The sixth 
lacuna of is that their Taylor expansion result (CI 3) grossly overestimates the convolution 
integral for large track lengths. 

An important observation must be made at this juncture. The above discussion along with 
Fig. 1 is the easiest way to check the inadequacy of the MT approximation (C13) based on Ai?, 
without getting bogged-down with experimental complications. If, however, actual RHIC data on 
jet quenching are to be fitted then due consideration must be given to the temperature T > 0.2 
GeV, plasma life time tuje < 10 fm/c, cylindrical geometry of Bjorken expansion, production 
configuration (r, 0) of the leading parton in the transverse plane, jet track lengths L(r, 0) < 10 fm 
etc. Even then, the detailed convolution integral should be computed using Simpson's quadrature 
at every stage rather than the AE' approximation. 



Ratio R versus E at various L 
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FIG. 1: The ratio R 
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^'^'^ plotted as a function of parton energy E at various track lengths L. 



Simp 



The curve corresponding to L = 10 fm has not been shown because R shoots to 170 for £^ = 5 GeV. 



Conclusions 



In view of the recent RHIC experiments the subject of jet quenching has become a very impor- 
tant research area because one must know the precise energy loss suffered by the leading parton 
due to the coUisional encounters and radiative bremsstrahlung. For calculating the coUisional 
energy loss one should start from a FP equation consistent with special relativity and work out 
all properties of the distribution function D(p,t) in conceptually valid manner. In this context 
refs. yj,|20 give excellent physical picture of the mechanisms but, unfortunately, most of their al- 
gebraic equations as well as numerical graphs seem to be unreliable. We will address these issues 
in detail in a subsequent paper. 
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